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We present evidence from computer simulation that the slowdown of relaxation of a standard 
Lennard-Jones glass-forming liquid and that of its reduction to a model with truncated pair poten- 
tials without attractive tails is quantitatively and qualitatively different in the viscous regime. The 
pair structure of the two models is however very similar. This finding, which appears to contra- 
dict the common view that the physics of dense liquids is dominated by the steep repulsive forces 
between atoms, is characterized in detail, and its consequences are explored. Beyond the role of 
attractive forces themselves, a key aspect in explaining the differences in the dynamical behavior of 
the two models is the truncation of the interaction potentials beyond a cutoff at typical interatomic 
distance. This leads us to question the ability of the jamming scenario to describe the physics of 
glass-forming liquids and polymers. 

PACS numbers: 05.10.-a, 05.20.Jj, 64.70.kj 



I. INTRODUCTION 

Since van der Waals, the statistical mechanics of liq- 
uids has been built on the idea that the structure of dense 
liquids is predominantly determined by the steep repul- 
sive forces between molecules and that longer-ranged at- 
tractive interactions merely establish a cohesive back- 
ground that affects the thermodynamics, but neither the 
structure nor the dynamics [l|-|4| • The physical rationale 
for this idea is that, in a dense liquid, density fluctuations 
coming from the displacement of the atoms or molecules 
involve the short-ranged forces. Longer-ranged forces, 
which are associated with a large number of particles, 
average to zero and have thus a negligible contribution to 
these fluctuations. This idea has been incorporated into 
an operational scheme by Weeks, Chandler and Ander- 
sen (WCA) [3| who proposed to split the intermolecular 
potential into a short-ranged repulsive component and a 
longer-ranged, more slowly varying component, and to 
treat the latter as a perturbation. This has provided 
a powerful tool to study liquids. Especially for simple 
atomic liquids, it has led to a very successful description 
of the structure, at the level of the static pair density 
correlations, and of the thermodynamics [1, 

As already understood in part by van der Waals,and 
stressed repeatedly in modern implementations 0, 0-01 > 
this picture which neglects the fluctuations associated 
with the attractive, and more generally long-ranged, in- 
teractions becomes less accurate or even breaks down in 
some cases. For instance, when the liquid has large-size 
inhomogeneities, such as an interface, or when the at- 
tractive forces are short-ranged and directional, the net 
vector force on a given molecule or on a liquid domain 
that results from attraction with the rest of the system 
no longer vanishes. Then, a homogeneous, mean-field ac- 
count of the attractive interactions is no longer valid. As 
is also well-known, on the low-density side of the liquid 



range, near the critical point or near the gas-liquid spin- 
odal induced by attractive forces, density fluctuations de- 
velop on all lengthscales and become large, which again 
makes the mean-field description inappropriate. 

If the van der Waals picture and the WCA theory of 
liquids have been extensively studied and tested as far 
as the pair structure and the thermodynamics are con- 
cerned, much less has been undertaken for the dynamics. 
Actually, the WCA theory and the associated division of 
the pair potential have not been submitted to systematic 
investigations [8, 9]. We have recently revisited this point 
in the context of glass formation and found that a binary 
Lennard-Jones model and its WCA reduction to the re- 
pulsive components of the pair potentials, when studied 
at the same density, show very similar static pair den- 
sity correlation functions but widely differing relaxation 
times as one lowers the temperature towards the glass 
phase [l(]| . In the present article, we expand on this first 
study and present a more exhaustive set of results for 
both the structure and the dynamics of the two liquid 
models. It confirms that a simple atomic glass-forming 
liquid in its viscous regime is manifestly an exception 
to the van der Waals picture for what concerns its dy- 
namics, but not its pair structure. Although we cannot 
unambiguously assign this observation to a unique physi- 
cal cause, it seems likely that it indicates that the viscous 
liquid is more heterogeneous than the mere consideration 
of its pair density correlation function would suggest. 

In the previous short account of this work fioj , we have 
described the effect of the attractive forces, included in 
the full Lennard-Jones model but not in its WCA reduc- 
tion, as "nonperturbative" in the viscous liquid range. 
We address here more thoroughly this issue. We try to 
characterize the qualitative difference in dynamical be- 
havior between the two models, by studying in particu- 
lar the scaling of the relaxation time with density and 
temperature. This scaling appears as a ubiquitous (but 
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approximate) property of glassforming liquids and poly- 
mers [UHlH]; it is obeyed by the Lennard- Jones model 
but is strongly violated in the purely repulsive WCA 
model. 

The presence or absence of density scaling in the dy- 
namics of the binary Lennard- Jones liquid and its WCA 
reduction is clearly the consequence of the presence or 
absence of an attractive tail in the pair potentials of the 
two models. However, one should be cautious about the 
conclusion to be drawn from this fact. It is actually not 
obvious whether this is due to the attractive nature of the 
tail per se or to the introduction of a cutoff at a typical 
interatomic distance in the pair potentials, irrespective of 
the attractive or repulsive nature of the neglected longer- 
ranged interactions. Density scaling is indeed an exact 
property of purely repulsive power-law potentials, and it 
has been recently shown that both the structure and the 
dynamics of the binary Lennard- Jones model can be re- 
produced by a binary mixture of atoms with purely repul- 
sive power-law interactions [TH, On the other hand, 
as we shall also stress, density scaling is strongly violated, 
at low temperature and pressure, in systems of repulsive 
harmonic-like spheres, for which the pair potentials van- 
ish beyond the typical interatomic distance [l?], EU • 

Building on these observations, we address the rel- 
evance of the "jamming scenario" 19] for describing 
the glass transition of liquids in their experimentally 
accessible range of density. This leads us to ques- 
tion the recently suggested equivalence between jamming 
phenomenon and colloidal glass transition on the one 
hand and glass transition of supercooled liquids on the 
other 0. 

In Sec. [IT] we describe the technical detail of our work. 
In Sec. lIIII we contrast the static and dynamic behavior of 
the two models under study. We study the scaling with 
density and temperature of the dynamics of the models 
in Sec. IIVI and we address the role of attractive forces 
per se versus that of introducing a cutoff at a typical 
interatomic separation in the pair potential in Sec. [V] 
Sec. IVII is devoted to a discussion of the connection be- 
tween glass-forming liquids and systems near jamming. 
Finally, we conclude the paper in Sec. IVII1 



II. MODELS, SIMULATION, AND PHASE 
DIAGRAM 



A. Models 

We compare the structure and the dynamics of a stan- 
dard three-dimensional model of glass-forming liquid, the 
Kob- Andersen 80:20 binary Lennard- Jones mixture [2lj |. 
and of its reduction to the purely repulsive part of the 
pair potentials proposed by WCA. In what follows, the 
former will be denoted by "LJ" and the latter by "WCA" . 
The interatomic pair potential between species a and /?, 
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FIG. 1: The P{T) equation of state for the LJ (filled symbols) 
and WCA (open sybols) models at different densities merely 
differ by a constant. 



with a, /3 — A, B, is given in the two systems by 



where 



4e 



a/3 



for r < r'i 



0, for r >r c a 



(1) 



a0 



is equal to the position of the minimum of 
v a p(r) for the WCA model and to a conventional cut- 
off of 2.5a a p (merely introduced for practical reasons 
with no impact on the physical quantities) for the stan- 
dard LJ model; C Q ^ is a constant that is fixed such that 
v aP ( r a/3) = 0- The difference between the LJ and WCA 
potentials is a purely attractive contribution, dominated 
by the r -6 contribution at large distance. 

We have performed Molecular Dynamics simulations 
in the NVE ensemble (after equilibration at a chosen 
temperature) with N — 900 — 1300 particles (depending 
on density) and we have studied a broad range of density 
p from 1.1 to 1.8. Lengths, temperatures and times are 
given in units of (Jaa, ^aa/^b, an d (ma\ A i '48e aa)' L ^ 2 re- 
spectively. It should be noted that the WCA truncated 
potential is continuous at the cutoff, but the resulting 
forces are not, which leads to a slow drift of the total en- 
ergy during very long simulation runs in the microcanon- 
ical ensemble. To cure this problem without introducing 
random collisions or a thermostat, we periodically rescale 
velocities with a very low frequency to maintain the total 
energy constant [22] ]. 

To check the robustness of our results, we have also 
considered the two-dimensional version of the L J and the 
WCA models, using a 65:35 mixture, because it is known 
to be less prone to crystallization in two dimensions than 
the 80:20 mixture Hi. 



B. Thermodynamics and phase diagram 

In the WCA theory, the repulsive truncated model is 
taken as a reference for the full LJ model, from which 
one can determine the structure of the liquid, whereas 



3 




FIG. 2: Phase diagram of the system with: the melting line of 
the monodisperse LJ system [24|], the coexistence curve [25l |. 
the spinodal line [2f| and the onset temperatures T onS et dis- 
cussed in Sec. IIVI of the binary LJ mixture, the T — ► jam- 
ming limit of the WCA model; the mode-coupling transition 
line T c (p) as fitted from the numerical data [271 ]. 



the attractive component is treated as a perturbation 
that allows one to compute thermodynamic quantities. 
Within this description, the structure and the dynamics 
of the LJ and WCA models should be compared at the 
same (p,T) state point. Their pressure then differs, with 
the attractive interaction roughly providing a negative 
background term. This is illustrated in Fig. [T] where we 
plot the simulation results for the equation of state, pres- 
sure versus temperature P(T), of the two models at four 
densities from p = 1.1 to p = 1.6: the pressure of the 
WCA model is roughly shifted up by a p-dependent con- 
stant from that of the full LJ model. As can be seen, the 
pressure of the full LJ mixture becomes negative at the 
lowest temperatures for p = 1.1, which is the sign that 
the liquid is in a metastable state in the two-phase re- 
gion inside the gas-liquid coexistence curve. The purely 
repulsive WCA model on the other hand has always a 
positive pressure and no coexistence region. 

To help characterizing the different state points that 
will be considered, we give in Fig. [2] the phase diagram 
in the (p, T) plane of the binary LJ model. For illus- 
trating the relevant range for the dynamics, we also plot 
on this diagram the estimated T c of the mode-coupling 
theory as obtained from a fit to the simulation data [13] ■ 
For the purpose of this work, the merit of such a fit is 
immaterial, the T c line being a convenient indicator of a 
temperature scale where dynamics has slowed down by 
about 4 decades as compared to the high-temperature 
fluid. The mode-coupling line ends up for p = 1.1 in 
the metastable liquid, as mentioned above, whereas den- 
sities p = 1.6 and p = 1.8 are clearly outside the conven- 
tional liquid range (one is always much above the criti- 
cal temperature). These densities actually correspond to 
conditions that, although easily investigated in Simula- 
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FIG. 3: Pair correlation functions of the (total) density fluc- 
tuations g(r) and of the concentration fluctuations g cc (r) for 
the LJ and WCA liquids at high and low temperature for 
the canonical liquid density p = 1.2. The different curves are 
vertically shifted, for clarity. 



tions of simple models, are not experimentally relevant, 
even in high-pressure experiments on glass-forming liq- 
uids (see also Ref. (28j). The density p = 1.2 is the 
standard condition at which the slowdown of relaxation 
of the binary Lennard- Jones mixture is usually consid- 
ered [2lj]. It corresponds to typical liquid (and viscous 
liquid) states. Most of our simulations are performed in 
the viscous regime in the region delimited by the onset 
temperature r onsot and the mode-coupling lines in Fig. [3] 
As for the WCA model, because of the absence of at- 
tractive interactions, it does not have a gas-liquid tran- 
sition and could in principle be studied at all densities 
down to zero temperature. In practice however, we found 
that the system crystallizes in simulations when cooled 
down to small but nonzero T at densities below p ~ 1.1, 
and this is the reason which prevents us from extending 
the mode-coupling line down to T = 0. On the other 
hand, directly at T = 0, the WCA model could in prin- 
ciple be subjected to various ".jamming" protocols with 
density as the control variable [2{|. We shall come back 
to these aspects in the following. 



III. CONTRASTING THE BEHAVIOR OF THE 
PAIR STRUCTURE AND OF THE DYNAMICS 

A. Pair correlation functions 

We first consider the structure of the two liquid mod- 
els, as characterized by the static two-body density cor- 
relations. We focus here on the pair correlation functions 
gaf}{ r )i the static structure factors S a p(q) having already 
been displayed elsewhere [jj . 

In Ref. [l(|, we have shown that the pair correlation 
function of the (total) density fluctuations, 

g(r) = x\ gAA(r) + 2 x A x B 9ab{t) + x 2 B ffBfl(r), (2) 
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where x a denotes the concentration of species a, is ex- 
tremely similar for the LJ and the WCA models and that 
its shape varies only weakly with temperature. This is 
illustrated in Fig. [3] for the canonical density p = 1.2. 
Weeks, Chandler and Andersen Q pointed out that, in 
the case of a binary mixture, the coincidence between 
the structure of the two models should not be as good 
for the pair correlations involving the fluctuations of con- 
centration (5c = xbSpa — xa^pb), as those are less con- 
strained at high density than the (total) density fluctua- 
tions (8p — xaSpa + xbSpb)- We have checked this by 
plotting the pair correlation function of the concentration 
fluctuations, g C c{r), defined as 

g cc (r) = x 2 B g A A{r) - 2 x A x B gAs{r) +x 2 A g B B(r). (3) 

The results are also shown in Fig. [3] As anticipated, the 
difference between the LJ and WCA models is still weak 
but more pronounced for g C c(r) than for g(r). 

Finally, we display in Fig. U the evolution with density 
of the partial pair correlation functions g a p{r) for the 
two liquid models. As there is less difference between the 
two models at high temperature, we concentrate on the 
low temperatures, thus emphasizing those state points 
where differences are more pronounced. At high density, 
all curves nearly superimpose, but one can see notable 
differences between the two models at p = 1.1, espe- 
cially for the correlation functions involving the minority 



species B. At such a density, attractive interactions start 
playing a nonnegligible role as the full LJ system enters 
the metastable liquid region inside the coexistence curve 
and approaches the spinodal where density fluctuations 
grow very large. 

It is perharps not surprising that BB correlations are 
more affected by removing the attractive part of the po- 
tential, as the interaction parameters of the model were 
specifically chosen to favor attraction between AB par- 
ticles, such that B particles can efficiently frustrate the 
crystallization of the majority A component. By remov- 
ing the attractive interactions, we see that the first peak 
of gBB^r) is more pronounced for the WCA model, while 
the first peak in gAB (?) is less pronounced, showing that 
the effect engineered by Kob and Andersen is less effi- 
cient. This also explains why the present WCA mixture 
crystallizes more easily at low density than its LJ coun- 
terpart. Note finally that gBB{r) dominates the behav- 
ior of gcc(f), see Eq. ((3]), which explains why g cc (r) is a 
more sensitive probe of the structural differences of the 
two models than <?(r), which is dominated by gAA (?")■ 

Although we discussed in detail these small differences, 
our results on the whole provide one more confirmation 
of the validity of the WCA theory for the equilibrium 
pair structure of dense liquids. The rather small differ- 
ences observed between the LJ and WCA models could 
probably be captured by a perturbative treatment of the 
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FIG. 5: Comparison of the time dependence of the self- 
intermediate scattering function F„(q,t) for qaAA — 7.2 of 
the LJ and WCA models at the same temperatures, shown 
for different densities. 



attractive tail, for instance along the lines suggested in 
Ref. 301. 



B. Dynamic correlations and timescales 

We now turn to the results concerning the dynamics of 
the two liquid models. As already shown in our previous 
publication [10j , the idea that the dynamical behavior of 
a liquid is essentially insensitive to longer-ranged attrac- 
tive interactions breaks down in an increasingly spectac- 
ular manner as temperature decreases. The difference 
between the liquid in the presence and in the absence of 
the attractive tail is small but already noticeable at high 




0.35 0.7 
l/T 



FIG. 6: Arrhenius plot of the relaxation time of the LJ (filled 
symbols) and WCA (open symbols) models for several den- 
sities. Note that the temperature range is different for each 
panel. 



temperature in the "normal" range (see also Refs. HQ) 
and it rapidly builds up in the low-temperature regime 
relevant for the glass transition phenomenon Here, 
we analyze in more detail the influence of the density. 

To analyze the dynamics of the models, we have mostly 
studied the time dependence of the self-intermediate scat- 
tering functions defined by 



F s (q,t) 



N \ 

y^ e *q-Mt)-ri(o)) 



(4) 



with qaAA — 7.2, which we kept constant for all state 
points. This corresponds roughly to studying single par- 
ticle displacements over a length scale comparable to the 
inter-particle distance. Representative data points are 
shown in Fig. [5] for three densities and a broad range 
of temperatures. These figures strikingly illustrates how 
the small differences in the relaxation dynamics of the 
high-tcmperature liquids become a dramatic effect at low 
temperatures. For p = 1.2 and T = 0.45 for instance, the 
LJ system hardly relaxes in the timescale of the simula- 
tion while this temperature corresponds to a modestly 
supercooled state. When p increases, these differences 
decrease, but note that even for p = 1.6 and T = 1.8 (a 
temperature which is 50 % larger than the critical point, 
see Fig. ^ the time correlation functions still differ by a 
large factor. 

From the time decay of the self-intermediate scatter- 
ing function, we extract a relaxation time, r, defined in 
practice as F s (q, t) = exp(— 1). We show in Fig.[S]an Ar- 
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FIG. 7: Two-dimensional 65:35 Kob- Andersen binary LJ mix- 
ture and its WCA reduction: Top: Partial pair correlation 
functions g a p(r) for a typical liquid density p — 1.2 at low 
temperature. Bottom: Time dependence of the self interme- 
diate scattering functions under the same conditions. 



rhenius plot of this relaxation time for several densities. 
The figure clearly illustrates that the dramatic difference 
between the two models that is seen at liquid densities 
(with already a difference of 3 orders of magnitude at the 
lowest temperature at which we can equilibrate the full 
binary LJ model) decreases slowly with density. In such 
a logarithmic representation covering many decades, the 
presence or absence of the attractive tail of the potentials 
appears irrelevant at very high densities only. However, 
as stressed above, such densities are not realistic for ac- 
tual glass-forming liquids. 



Jones forces has a dramatic impact on the dynamics can 
in fact be drawn by comparing the results already pub- 
lished in the literature for different Lennard- Jones mix- 
tures [U [3l[ and for their WCA truncations (32l, l33j . 

Additionnally, since their initial publication 10], our 
results for the three-dimensional Kob- Andersen 80:20 
mixture have now been confirmed by independent stud- 
ies El El. 



D. To be, or not to be (perturbative) 

Before closing this section, we would like to discuss 
how to best characterize the difference in the dynamical 
behavior of the LJ and WCA models at liquid densities 
and what to make of it. In our earlier publication [To| . 
we have described the effect of the attractive tail of the 
pair potentials on the dynamics as "nonperturbative" . At 
a quantitative level, the gap between the characteristic 
time scales of the two models indeed becomes enormous 
and widens as temperature decreases. By itself, this ob- 
servation undoubtedly invalids any claim that attractive 
tails can be neglected when describing the dynamics (the 
relevance, or not, of the attractive nature of the tail will 
be discussed in section [VJ, but it is not enough to justify 
the qualifier "nonperturbative". In principle, this quan- 
titative difference, no matter how large, could be scaled 
out through the introduction of parameters which them- 
selves could be computed via a perturbative treatment 
of the longer-ranged interactions. This would be true for 
instance if the dynamics could be predicted on the sole 
basis of the knowledge of the pair structure, with mi- 
nor effects in the latter being strongl y a mplified in the 
former. However, we show elsewhere [34| that such ap- 
proaches are unable to account for the diverging behavior 
of the LJ and WCA liquids in the viscous regime as tem- 
perature is lowered. In addition, we have found evidence 
that the difference in the dynamics of the two systems 
is not only quantitative but also qualitative. This is the 
point we now discuss in more detail. 



IV. ON THE DENSITY SCALING FOR THE 
RELAXATION TIME 



C. Two-dimensional case and other examples 



Density scaling in supercooled liquids and 
polymers 



We note in passing that the pattern that we have found 
here, namely that the WCA reduction of the binary LJ 
mixture to a truncated purely repulsive systems has vir- 
tually no effect on the pair correlations of the total den- 
sity but strongly affects the dynamics, is not unique. 
We have also obtained the same effect for instance in a 
two-dimensional version of the binary 65:35 LJ mixture 
model, as illustrated in Fig. [7] 

We point out that, although not appreciated before, 
the conclusion that the WCA truncation of Lennard- 



The simplest way for trying to scale out the difference 
in the T-dependences of the LJ and WCA liquids is to 
renormalize the temperature by one density-dependent 
parameter chosen to make the curves of the two mod- 
els coincide at high temperature, where a perturbative 
treatment of the attractive tails more likely to apply. 

The most convenient way to proceed is to fit the high-T 
data to an Arrhenius form, 



r(p.T) 



3 (p)exp[ J B 00 (p)/T], 



(5) 
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which indeed provides a good empirical description, ir- 
respective of the physical meaning one is willing to put 
in its use. This procedure suggests to rescale T for the 
whole range under study by the extracted energy param- 
eter Eggip). When applied to glass-forming liquids and 
polymers [11], this procedure has been shown to provide 
a very good collapse of the relaxation and viscosity data 
at different densities, according to 



EUp) 

T 



(6) 



With specific functional forms for the density-dependent 
scaling parameter (see also below) , this relation has been 
successfully applied to an impressive variety of glass- 
forming liquids and polymers [llTflil ]. Thus, there is 
enough empirical evidence to take this scaling property 
as a genuine characteristic of glass formation in these 
systems. 

It must be noted that Eq. ([6]) is far from being triv- 
ial, as it connects the dynamics in the high-temperature 
liquids to the behaviour found in the viscous, low tem- 
perature regime. More precisely, Eq. ^ shows that two 
physical quantities are in fact proportional: the high- 
temperature activation energy, Eoo , and the "onset tem- 
perature" for slow dynamics, T onsot , which could be oper- 
ationally defined at the point where T\{x) departs from 
its high-T (Arrhenius) behavior: 



Tonset(p) OC E oc (p). 



(7) 



Physically, T onse t marks the temperature below which the 
liquid starts behaving in a more collective and heteroge- 
neous manner [35l - l38| . 

As expected, and already displayed in Ref. [Toj ] . the 
data for the full LJ binary mixture obey very well the 
density scaling of Eq. © , see Fig. 2a of Ref. [10] (a small 
deviation from the data collapse can be observed at the 
lowest density p = 1.1 corresponding to the metastablc 
liquid inside the gas-liquid coexistence curve). On the 
other hand, the scaling relation is strongly violated for 
the WCA liquid: see Fig. 2b of Ref. [Io|. 

One important consequence is that the isochoric 
fragility, which characterizes the steepness of the T- 
driven slowdown of relaxation and can be quantified by 
the derivative d log[r(p, T)/too] jd log T taken at constant 
p and evaluated at the glass transition (or for any given 
value of t/tqo), is essentially independent of density for 
the LJ model, as also found in real glass-forming liquids 
and polymers fill ], but strongly depends on density for 
the WCA model. This presence or absence of (even ap- 
proximate) density scaling represents, we claim, a quali- 
tative difference between the two liquid models. 



B. Empirical data collapse 

If one wishes to also collapse the relaxation data at 
all densities for the WCA model, one must introduce at 
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FIG. 8: Empirical collapse of the relaxation-time data for the 
LJ and WCA models at densities from p = 1.1 to p = 1.8. 
Top: Using Eq. ([9]) with parameters listed in Tabled] Bottom: 
Using Eq. Q with parameters listed in Table [I] The absolute 
scale of y(p) is chosen so that the deviation from the high-T 
behavior takes place around T>nsct = yEoo ~ 1. 



least one additional, p-dependent parameter. This pro- 
cedure may have only little physical significance per se, 
but it allows one to quantify and discuss more precisely 
the deviation from density scaling, which is obeyed by 
the LJ model only. 

We have found that the two functional forms, 



r(p,T) 



exp(£; oo (p)/T) 



v(p) 



eUp) 



T 



and 



r(p,T) 



x{p) 



EM 



T 



(8) 



(9) 



both lead to an excellent collapse on a master-curve of the 
L J and WCA data, as seen in Fig. [8] It is straightforward 
to show that the density scaling of Eq. ([6]) is recovered 
when either x{p) or y(p) is independent of density. 

We arrived at the formula in Eqs. © O by using the 
possibility that Eq. (J7J) is not valid and that one needs 
two distinct energy scales to describe the behaviour of 
the WCA model. This additional freedom can be incor- 
porated in a variety of ways, and we present two of them 
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Tonsct 




y 




1.1 


0.197 


0.885 


4.49 


0.586 


0.517 


1.6 


3.09 


1 


1.2 


0.513 


2. 


3.89 


0.685 


0.833 


2.55 


3.06 


1 


1.4 


1.51 


5.25 


3.47 


0.897 


1.75 


5.6 


3.20 


1 


1.6 


3.2 


10.6 


3.31 


0.958 


3.33 


10.9 


3.27 


1 


1.8 


5.6 


18. 


3.21 


1.0 


5.67 


18 


3.17 


1 



TABLE I: The parameters for the empirical data collapse. 
The parameter y is equal to the ratio between Boo and Tonset- 
It fluctuates around w 3.2 for LJ at all densities, and decreases 
strongly with p for the WCA. 



in Fig. [8] In a recent work 39] , our published data [HI 
for the LJ and WCA models have been re-analyzed and 
collapsed for all densities by means of a specific functional 
form which similarly makes use of two independent en- 
ergy scales; however, this collapse only holds in the vis- 
cous regime, while we also consider the high-temperature 
regime in Fig. [3] 

As expected, we find that for the LJ model, the two ad- 
ditional parameters y(p) and x{p) are indeed essentially 
independent of density. While x = 1 is used in Fig. [51 
the y values are given in Table U with small variations 
from a constant that have no systematic dependence on 
density within statistical accuracy. On the contrary, to 
collapse also the data for the WCA model, significant de- 
viations from constant behavior (with systematic trends 
as a function of density) are needed for both x(p) and 
y(p), see Table HI 

The presence or absence of density scaling for LJ and 
WCA models is a qualitative effect, which we have em- 
pirically quantified by introducing the ^parameter x or y 
in Eqs. ((HI EI)- As pointed out in Ref. [3E|, the differences 
found in the parameters obtained for the two models are 
modest (x or y vary by about 50 %) compared to the 
observed differences in the relaxation times themselves; 
but this is of course expected for parameters characteriz- 
ing the temperature dependence of the logarithm of the 
relaxation times, and this does not contradict our con- 
clusion that the presence or absence of density scaling 
represents a qualitative feature of the slowing down of a 
glass-former. 

From Eq. (j8j we see that a sensible definition of the 
onset temperature becomes 



Tonset(p) OC y(p)E oa (p). 



(10) 



This is nicely illustrated in Fig. [5] (bottom) where T onset 
can be estimated as the point at which the relaxation 
data starts to significantly deviate from the high-T be- 
havior; the resulting values of T onsct are given in Ta- 
ble [T] The physical significance of the values obtained in 
this way are independently confirmed by checking that 
Tonsct also corresponds to the crossover at which the 
T-dependences of the self-diffusion constant, D sc u{p,T), 
and of the relaxation data, r(p,T), begin to split, see 
Fig. IH1 for an illustration at p = 1.2. These values are 
independently confirmed in Ref. [40]. 



1 0.04 




FIG. 9: Decoupling between the T dependences of the self dif- 
fusion constant -D so if and of the relaxation time r for p = 1.2 
as a function of the temperature T rescaled by the indepen- 
dently measured T onse t . For both models, the product D 3e ll t 
begins to increase below T « T onsc t, which differ by a factor 
1.6 for both models, see Table [I] 



The data compiled in Table fl] actually imply that the 
LJ and WCA models have distinct dynamics in the high- 
temperature regime already since of both models 
only become equal at very high density. Both models are 
also characterized by distinct onset temperatures, which 
differ by a factor « 1.6 at p = 1.2, and, again, slowly con- 
verge to similar values at large density. Finally, the qual- 
itative character of the impact of the truncation of the 
potential on the slow dynamics is illustrated by the fact 
that the isochoric fragility of the L J model is independent 
of density while fragility strongly increases with density 
for the WCA model, as parametrized by the evolution 
of either x{p) or y(p) which is qualitatively different for 
both models. 



C. Analogy with power-law repulsive spheres 

In the case of the LJ liquid, the density scaling can be 
expressed, with essentially as good data collapse as with 
Eq. ([6]) , by using a power-law density dependence for the 
scaling parameter (see Fig. 2c of Ref. [10| showing that 
Eoo(p) ~ P 5 ), namely, 



T 



(11) 



with 7 ~ 5. In this form, density scaling was argued to be 
related to the property that LJ models are "strongly cor- 
relating liquids" 41], which means that the fluctuations 
in the potential energy, 



U = ^2v(\Ti- rj \), 



1<J 



and those in the virial part of the pressure, 
W = -^5>(|r,- rj -|), 



(12) 



(13) 
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with w(r) — rv'(r), are correlated at all state points. The 
correlation is characterized by a single parameter, T, de- 
fined by a linear fit through a scatter plot of the time fluc- 
tuations of U and W. The connection between the den- 
sity scaling of the dynamics and the strong-correlation 
property of the fluctuations shows up in the approximate 
relation: 7 ~ T. 

The rationale which has first been put forward to jus- 
tify both properties 0, [HI 5H is as follows. A dense LJ 
liquid is dominated by the repulsive forces (the van der 
Waals picture of liquids), but the relevant part of the re- 
pulsive interaction involves interatomic distances in the 
range near to (but less than) the minimum of the pair 
potential; in this range, the potential looks like an effec- 
tive power law r -37 with an effective repulsion which is 
steeper than the bare r~ 12 repulsive component of the 
LJ potential, i.e. 37 > 12. In this picture, the LJ model 
essentially behaves as a repulsive r -37 power-law inter- 
acting model, a model that is known to rigorously lead to 
the two considered properties, with indeed a strict equal- 
ity 7 = r. 

The simple rationale given above has however a major 
flaw: it predicts that the WCA and the full LJ mod- 
els should behave in exactly the same way as their pair 
potentials coincide for distances shorter than the min- 
imum. This, as shown above, is clearly wrong. Thus, 
the WCA-like argument that the fluctuations in dense 
Lennard- Jones liquids (density scaling and U — W cor- 
relation) can be understood on the basis of the steep 
repulsive core only is, in general, not valid. 

An improvement over the argument was later sug- 
gested by Dyre and coworkers [421. They argued that 
the relevant interatomic distances actually include the 
vicinity of the minimum on both sides of it and that a 
better description of the LJ potential is through a power 
law complemented by a linear term; as the latter has es- 
sentially no effect on the fluctuations at constant density, 
the fluctuation properties of the full LJ system are thus 
predicted to be those of a properly adjusted power-law 
repulsive potential. The WCA truncated potential being 
indeed less well described by a power law plus a linear 
term (for distances above the minimum), this argument 
could allow one to sidestep the criticism raised above. 

We have repeated the analysis of the correlation be- 
tween the fluctuations of the potential energy and those 
of the virial for the two models. We find that, as expected 
from the work of Dyre and coworkers, the LJ model has a 
strong correlation at any given density. The WCA model 
also has strong correlations between energy and virial at 
any given (p,T) state point [43j, as illustrated in Fig. [TU1 
for p = 1.2 and T = 0.45. 

The parameter F that we extract from the slope of the 
correlation plot is reported in Fig. [TUJ for a broad range 
of state points. It is virtually independent of temper- 
ature and weakly dependent on density in the LJ case, 
but has a stronger density dependence (and temperature 
dependence when p = 1.2) in the WCA case (see also 
Ref. [Hj]). Therefore, if one neglects the modest density 
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FIG. 10: Top: Linear correlation between the fluctuations of 
the virial W and the fluctuations of the potential energy U 
for the WCA liquid at p = 1.2 and T = 0.45, with a slope 
F. Bottom: Variation with temperature of the slope F of the 
correlation between the fluctuations of W and U for the LJ 
and the WCA liquids at several densities between 1.2 and 1.8. 



dependence of T for the LJ model (in Fig.QJil the effective 
slope r varies between 4.5 and 5.2 over a wide range of 
density and temperature), one can indeed conclude that 
the latter is "strongly correlating" over the whole liquid 
range with an effective, constant, value of T of about 5. 
Of course, the data in Fig. [TU] underlie that there is still 
some arbitrariness in the choice of the value of T, and 
the quality of the density scaling depends on this choice. 

On the other hand, as shown in Fig.QJJJ the variation of 
T in the WCA liquid is probably too large to be replaced 
by a single, effective constant. In this case, neither the 
fluctuations nor the dynamics are purely determined by 
reducing the pair potentials to repulsive power-law ones 
with an exponent chosen to reproduce the effective steep- 
ness of the potential near the minimum. Thus, although 
both phenomena (density scaling and strong-correlation 
property) seem connected to each other, the example of 
the two models presented in this study shows that it is 
not obvious to know a priori for which materials these 
observations may apply or not. 
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V. IS THE ABSENCE OR PRESENCE OF 
(APPROXIMATE) SCALING DUE TO 
ATTRACTION OR TO TRUNCATION IN THE 
POTENTIAL ? 

We have seen that the presence (in the LJ model) or the 
absence (in the WCA model) of an attractive tail in the 
pair potentials has a large quantitative influence on the 
dynamics and strong consequences for the fluctuations 
and the relaxation properties of the liquids. In partic- 
ular the (approximate) density scaling observed for the 
temperature dependence of the relaxation time in glass- 
forming liquids and polymers (with the resulting density 
independence of the isochoric fragility) is found in the 
LJ model but not in the WCA one. A central question 
then to be raised is whether the differences stem from the 
attractive character of the tail per se or from the trunca- 
tion of the range of the potential to typical interatomic 
distances. 

A first hint that truncation is the key feature is pro- 
vided by looking at the behavior of systems of spheres 
with purely repulsive power-law potentials. As already 
mentioned, such liquids are (rigorously) "strongly corre- 
lating" and show an exact density scaling of the relax- 
ation time, as a consequence of the scale-free power law 
behavior of the potential. This example evidences that 
the presence of attractive interactions is not a necessary 
ingredient in establishing these properties. 

Additional evidence along the same lines has recently 
been provided by Pedersen et al. (l6| who showed that 
the pair structure and the dynamics of the binary LJ 
mixture at liquid densities can be very well reproduced by 
replacing the LJ pair potentials by nontruncated power- 
law repulsive potentials with an appropriately adjusted 
exponent. Note that they chose a repulsion of the form 
r -i5.48 guch that 7 = r = 5.16, in rough agreement with 
the data shown in previous sections. 

To confirm that truncating the potentials beyond a 
cutoff of the order of typical interatomic distances is re- 
sponsible for the absence of density scaling seen in the 
WCA model, we revisit the example of fluid mixtures of 
repulsive harmonic spheres, with pair potentials 

v a p(r) = — — ) ,forr<cr Q/3 

s V <*ap) (14) 

= 0, for r > <T a/ 3, 

where s — 2; these systems are commonly used in the 
context of zero-temperature jamming phenomena [2^ |. 
A wide range of temperature and density was studied by 
computer simulation in Refs. [UlSl- As noticed in these 
references, the isochoric fragility of the system strongly 
depends on density. We have replotted in Fig. [TT] the 
relaxation-time data in the way described in section IV 
(and used in Fig. 2 of Ref. [l(| ) by rescaling the temper- 
ature with a density-dependent energy parameter E x (p) 
chosen to make the high-T data collapse on a single curve. 
This plot clearly shows that, just like the WCA model 
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FIG. 11: Rescaling of the relaxation time data for the har- 
monic repulsive spheres at low temperatures after having nor- 
malized the temperature by the energy parameter E x {p) cho- 
sen to make all curves collapse at high-T. The density p is 
indicated in the figure. This evidences a dramatic breakdown 
of density scaling in harmonic spheres, similar to the one ob- 
served for the WCA model. 

and at odds with the full LJ one, there is no density 
scaling of the relaxation for truncated harmonic spheres. 

The similarity of behavior between the truncated WCA 
model and the truncated harmonic potential on the one 
hand and between the LJ model and the repulsive power- 
law model on the other is a strong indication that the 
presence or absence of density scaling in the dynamics 
of a glass-former results from the presence or absence of 
a truncation of all pair interactions at a typical inter- 
atomic distance. This observation would perhaps not be 
very surprising at much lower densities and temperatures 
(see Fig. [2]) since, by construction, the WCA potential co- 
incides with the one of harmonic spheres in Eq. (|14[) near 
the cutoff. The surprising feature is that this analogy 
seems to be relevant up to the liquid densities studied 
in the present work, and it results in the WCA model 
behaving differently from the LJ system. As discussed 
in the following section, this conclusion casts doubts on 
a description of supercooled liquids in terms of the jam- 
ming scenario. 



VI. GLASS TRANSITION VERSUS JAMMING 
PHENOMENON 

The jamming paradigm (l9l l44j has been put forward 
to bring together in a common picture a wide breadth 
of phenomena and systems involving sluggish dynamics 
and freezing in an amorphous state. A step to go be- 
yond qualitative comparisons has been taken with the 
proposal that the slowing down of all jamming systems, 
whether driven by temperature, density or applied force, 
is controlled by a zero temperature and zero applied force 
critical point, "point J" [29]. This proposal has been crit- 
icized on several grounds, in particular concerning the 
uniqueness of point J itself [45|, H(| , and the possibility 
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that a distinct glass critical point, which could be dubbed 
"point G" , controls the finite temperature dynamics of 
jamming systems [53] . Here, we focus on the ability of the 
jamming scenario and of the associated point J (whether 
unique or not) to describe the slowdown of relaxation of 
actual glass- forming liquids (and polymers). 

A transition near point J only takes place at zero tem- 
perature and moderate density, and it can only be rigor- 
ously defined for interaction potentials that are zero be- 
yond some distance and repulsive within it. Therefore, to 
be applicable to liquids in their (experimentally) acces- 
sible range, the jamming scenario must assume that the 
role of the attractive tail and more generally of longer- 
ranged (longer than the typical interatomic distance) in- 
teractions is negligible in the physics of the slowing down. 
To extrapolate the behavior of a real liquid to such con- 
ditions, one must then get rid of the attractive forces, 
which give rise to the gas-liquid coexistence curve and 
prevent one from taking the liquid to zero temperature 
below some density (see the phase diagram in Fig. [2]), 
and eliminate the pair interactions beyond some cutoff 
that physically determines the typical inter-particle dis- 
tance. Note that this is precisely what the WCA pro- 
cedure achieves. However, as we have seen above, the 
behavior of the WCA model strongly deviates, both qual- 
itatively and quantitatively, from that of the LJ model 
and more generally that of realistic liquids, except at 
very high densities that are beyond any experimentally 
accessible range. It is this observation which leads to 
question the ability of the jamming scenario to describe 
the dynamics of viscous liquids. 

Actually, systems of particles interacting through trun- 
cated repulsive interactions, such as the harmonic-like 
models in Eq. (| 14|) or the WCA model, behave as "ef- 
fective" hard-sphere systems at low temperature and low 
pressure. In such thermodynamic conditions, it is the 
behavior of the pair potential near the cutoff which gov- 
erns the physics, and if temperature is low, soft repulsive 
spheres can be seen as "disguised" hard spheres [13, lla - 
As shown in Fig. [TT] and discussed in detail in Ref. [17| 
this directly implies that the isochoric fragility of such 
liquids strongly depends on density. This is a first quali- 
tative distinction between the physics near the jamming 
transition and the one of supercooled liquids. 

A second consequence, discussed in Ref. [2(|, is that 
relaxation data for soft spheres with truncated repulsive 
interactions at different state points converge in the low 
T and low P limit to a hard-sphere-like scaling, such 
that [U, 



r(p,T) 



= H 



P(P,T) 



T 



(15) 



where Too has a residual (trivial) dependence on either T 
or P {e.g. in l/y/T or vP) that it is now important to 
take into account as the range of pressure and tempera- 
ture spanned is very large [compared to the liquid range 
considered for Eqs. (j MTTj) ]. The domain of application of 
this scaling is demonstrated in Fig. [T2] where we include 
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FIG. 12: Hard-sphere-like scaling of the relaxation time. We 
use Eq. (|15p and follow the evolution of the relaxation dynam- 
ics from hard spheres [48j |. to harmonic spheres [l7t | . to WCA 
and finally LJ systems, which exhibit increasingly larger de- 
viations from hard-sphere behavior. 



data obtained from simulations of hard spheres taken 
from Ref. (48|. The harmonic-sphere data of Fig. [TT] in- 
deed converge to the hard-sphere behavior when de nsity 
is not too large, in agreement with the results of Ref. [20j ]. 
and display small deviations from Eq. (1151) when density 
is larger. These deviations were empirically addressed in 
Ref. [13, [H| using a (p, T)-dependent "effective" volume 
fraction. Dynamics in such systems with truncated repul- 
sive potentials at low temperature and low-to-moderate 
pressure seems then dominated by (renormalized) free- 
volume or congestion effects, as for a hard-sphere fluid. 
Whether (or not) the slowdown of relaxation is controlled 
by a T = point J [29] or by a T > point G 0, S3, 
the jamming and glass transitions appear in any case 
strongly intertwined in this case. 

However, when moving from harmonic to WCA poten- 
tials in Fig. ll2[ we observe much stronger deviations from 
the hard-sphere scaling, and the data in fact converge to 
the vertical axis T/P = when density increases. This 
is expected, as in the large density limit we expect the 
relaxation-time data to be driven by the scaled variable 
T 5 / 4 /P, which becomes exact for the power law r re- 
pulsion. 

Making the final step from WCA to LJ models in 
Fig. [T^] we realize that the effect of the attractive forces 
becomes dramatic for normal liquid densities. Since at- 
tractive forces induce the presence of the gas-liquid phase 
separation, the pressure near the coexistence curve drops 
to zero when temperature decreases much more rapidly 
than in the absence of the attractive tails. The effect 
is indeed spectacular for the density p = 1.2 in Fig. Q2] 
since the LJ data have a qualitative behavior opposite to 
that of truncated purely repulsive systems. This echoes 
Voigtmann's remarks (28| emphasizing that WCA and 
LJ models appear more similar when using (p, T) than 
(P, T) variables, for the same reason. This observation 
is experimentally relevant as atmospheric-pressure data 
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for viscous liquids virtually all fall in this regime [28j, 
implying that the hard-sphere-like scaling in Eq. (fT5|) is 
not obeyed for liquids at normal conditions. 

This second argument confirms our claim that in the 
range of pressure, density, and temperature that corre- 
sponds to the actual liquid/supercooled liquid range, the 
glass transition of the liquid and that of jamming models 
described by truncated repulsive potentials are different 
and cannot be controlled by the same T — point J 
critical point. 

We finally make the aside that the models defined 
through truncated soft repulsive potentials differ when 
considered at high temperature and density. The (un- 
bounded) WCA potential leads to a system that con- 
verges to the full LJ model, as shown above; asymptot- 
ically, at large p and T, the WCA and the LJ models 
behave as a rescaled hard-sphere system (the rescaling 
is however not the same as at low pressure and temper- 
ature). On the contrary, the (bounded) harmonic-like 
models rapidly reach very high densities at which large 
numbers of atoms overlap, which seems to confer to these 
models some sort of mean-field behavior [49j , with pecu- 
liar physical consequences [50j . 

VII. CONCLUSION 

Through an extensive comparison of the behavior of 
a standard Lennard- Jones glass-forming liquid and that 
of its WCA reduction to a model with truncated pair 
potentials without attractive tails, we have shown that 
the slowdown of relaxation is quantitatively and qualita- 
tively different in the two models while the equilibrium 
pair structure remains very similar. The differences in 
the dynamics of the two models decrease as density in- 
creases but one has to reach unphysically large values 
of the density for seeing a full convergence of behavior. 
Clearly, the presence or absence of the attractive tails 
cannot be neglected in the viscous regime, where a liquid 
can no longer be simply described according to the con- 
ventional van der Waals picture: fluctuations associated 
with the tails of the interaction potentials do play a role. 

At this point, it should be stressed that, in spite of the 
observed differences, the purely repulsive WCA model is 
a glass-former, no less than the full LJ one. The attrac- 
tive tail is not a necessary ingredient to trigger a slow- 
down of relaxation. The WCA model displays features 
generically associated with glass formation. For example, 
as seen from Fig.[5J the associated time dependence of the 
self intermediate scattering function is characterized by a 
stretching and the appearance of a plateau as one lowers 
the temperature. Moreover, the dynamics is increasingly 
heterogeneous [3] and shows a decoupling between diffu- 
sion and structural relaxation, see Fig. |H1 Not surpris- 
ingly, the temperature dependence of the relaxation time 
data can also be fitted by the same functional forms, e.g. 



the Vogel-Fulcher-Tammann formula or a Bassler-type 
expression 39], for both models. This is enough to in- 
voke analogies, but not necessarily to assign the slowing 
down to the same mechanisms. 

The nonperturbative role of attractive forces may be 
tentatively attributed to the growth, as temperature de- 
creases, of some form of heterogeneities in the liquid that 
are not captured by the static pair density correlations. 
We have not investigated subtler characterizations of the 
liquid structure, e.g. via higher-order correlation func- 
tions or topological measures, but preliminary studies 
seem to indicate that these are more sensitive to the at- 
tractive tails than the pair density correlations and show 
structural differences between LJ and WCA models that 
grow as temperature is lowered [||], HH, [52j ■ This is cer- 
tainly a worthwhile line of research to pursue. 

In trying to characterize the physical significance of 
the increasing quantitative difference in the dynamics of 
the full LJ model and of the truncated repulsive WCA 
one, we have stressed that this difference could not be 
accounted for by a mere one-parameter rescaling of the 
data. (The fact that the difference cannot be reproduced 
either by a large amplification in the dynamics of the 
small differences in the static pair density correlations is 
considered elsewhere [13,|33]-) The temperature-driven 
slowdown of relaxation in the two models actually show 
qualitatively distinct behavior, as manifested by the ab- 
sence in the WCA model of the density scaling that is 
found in the full LJ model and in experimentally studied 
glass-forming liquids and polymers. Moreover, in the LJ 
and WCA models, this feature appears to be related to 
the properties of the correlation between the fluctuations 
of the virial and of the potential energy. 

Still concerning the physical significance of our results, 
we have addressed the question of whether the observed 
differences between the WCA and the LJ models in the 
viscous regime are due to the attractive nature of the 
missing tail (in the WCA model) or to the fact that the 
tail is truncated at a cutoff corresponding to typical inter- 
atomic distances (in the WCA model). We have provided 
strong evidence that the key feature is the truncation 
of the interaction potentials and we have stressed that 
models with truncated repulsive potentials show quali- 
tatively distinct temperature-driven slowing down than 
glass-forming liquids. This casts serious doubts on the 
possibility of describing glass formation in liquids by a 
jamming scenario in which the slowdown of relaxation 
is controlled by a critical point at zero temperature and 
density less than liquid densities. 
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